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ANALYTICITY OF THE DENSITY OF ELECTRONIC 

WAVEFUN CTIONS 

S. FOURNAIS, M. HOFFMANN-OSTENHOF, T. HOFFMANN-OSTENHOF, 
AND T. 0STERGAARD S0RENSEN 


Abstract. We prove that the electronic densities of atomic and 
molecular eigenfunctions are real analytic in R 3 away from the 
nuclei. 


1. Introduction and statement of the results 

We consider an TV-electron molecule with L fixed nuclei. The non- 
relativistic Hamiltonian of the molecule is given by 


H = H N)L ( R,Z) 


N 


£ 



+ £ 

l<i<j<7V 



£ 

1 <l<k<L 


Z\Z k 

Ri — Rk 


(l.l) 


where R = (Ri, R 2 , ■ ■ ■, Rl) G R 3i , Ri 0 Rk for k ^ l, denote the 
positions of the L nuclei whose positive charges are given by Z = 
(Zi, Z 2 ,..., Z L ). The positions of the N electrons are denoted by 
x = (#i,X 2 , • • • ,xn) G M 3JV , where x 3 denotes the position of the j'tli 
electron in M 3 . For shortness, we will sometimes write 


H = — A + V(x), (1.2) 

where A = A j, is the 3A-dimensional Laplacian, and V is the 

complete (many-body) potential. It is a standard fact that H with 
domain 1F 2,2 (1R 3JV ) is selfadjoint. 

We consider eigenfunctions "0 of //, i.e. solutions 0 G L 2 (M 3Ar ) to 
the equation 

i/0 = £0, (1.3) 

with E G R. Since we describe electronic wave functions, and the elec¬ 
trons are Fermions, 0 has to transform according to certain irreducible 
representations of the symmetric group & N . However, our results are 
independent of this condition and we do not impose it. 
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Analyzing the spectrum of H and calculating (usually by some ap¬ 
proximation scheme) the eigenvalues E and the corresponding eigen¬ 
function^) 0 is the central theme of most of the investigations done 
by quantum chemists and physicists. For the interpretation of these 
investigations the eigenfunction 0 is much too complex being a func¬ 
tion of 31V variables and hence the one-electron density p(x) plays 
a prominent role. It is defined by 



| 0 (%) | 


(1.4) 


where we use the notation x, = (aq, ..., Xj-i, x, Xj + \,..., x^) and 
die,- = dx i... dxj-±dxj + i... dxN- 

The mathematical analysis of H has mainly centered around the op¬ 
erator theoretical point of view, see for instance 0, |J and references 
therein. The fact that ( |1.3| ) is an elliptic partial differential equation 
has not been exploited in such depth; so many questions which are nat¬ 
ural from a PDE point of view are not really understood. In particular 
regularity questions concerning 0 and p are natural and interesting. 
Note first that V is singular in 


E = < a; E 


d3 N 


N L 

nn 

i=1 1=1 


Xi 


N 

n 

1 <i<j 


X; 


X. 


l)=o 


(1.5) 


and real analytic in R 3JV \ E. Hence by standard methods of elliptic 
PDE, see for instance |5|, 0 is real analytic in R 3JV \ E. The first 
results concerning the regularity of 0 on all of R 3JV are due to Kato 
[6J. He showed that 0 is Lipschitz continuous and first formulated 
the well known cusp conditions, which describe the behaviour of an 
eigenfunction near the points where two particles are close to each 
other [|6j. Theorems II and lib]. See also the important paper by Simon 
H0| in which the Coulombic many particle potential V is identified as a 
special member of the so called Kato class and some results concerning 
the regularity of solutions of Schrodinger equations (equations of type 
(|1.3j ) with general V) are given. 

Regularity results concerning the Coulombic case extending Kato’s 
result have been obtained more recently in [£ij and 0; see also the more 
complete references therein to other results concerning regularity. 

There are now two related problems: 


i) Describe in more detail how the specific structure of the singu¬ 
larities in E turn up in the nonanalyticity of 0. Partial results 
can be found in the references cited above. 
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ii) Analyse the regularity properties of the one-electron density, 
defined in (|1.4|) , an object which has an immediate physical 
interpretation (see any textbook on quantum mechanics, for in¬ 
stance [Q) and enters all approximation schemes in a crucial 
way (Hartree-Fock, Thomas-Fermi, Density Functional Theo¬ 
ries etc). 


For the regularity questions concerning p defined in (|hj) it suffices to 
consider the (non-symmetrized) density p defined by 


p(x) = 


/R3JV-3 


(x, X 2 , • • • , Xn)\ 2 dx 2 • • • dx 


N • 


( 1 . 6 ) 


It is not clear a priori that p is real analytic away from the nuclei 
since in ( fOl ) one integrates over subsets of S where ip is not analytic. 
In two recent papers (]I| and [§]) the present authors have shown that 
p is smooth away from the positions of the nuclei (or in the case of an 
atom, away from the origin). The natural question is now whether p is 
real analytic away from the nuclei. This will be answered affirmatively 
in this paper. Of course in the proof of this result new difficulties arise, 
in particular all the estimates have to be much more explicit. 


Theorem 1.1. Let ip G L 2 (R 3N ) satisfy the equation 


Hip = Eip, 


with £e! and H given by (|1.1|). Let the density p be defined as in 

0 ). 

Then p is a real analytic function in R 3 \ {l?i,..., Rl}- 


Remark 1.2 (Atoms vs. molecules). In order to keep notation simple, 
we will only give the proof of Theorem |1.1| in the case of an atom. In 
this case we only have one nucleus, which we place at the origin, so the 
potential V is given by 


V = 



(1.7) 


The necessary modifications for the molecular case were indicated in the 
proof of the smoothness results in PJ. In the present proof of analyticity 
one has to make similar changes when working with molecules. 


Remark 1.3 (Density matrices). We get analogous results for the one 
electron density matrix 7 1 (x,x / ) and the 2-electron density p 2 (x,x'), 
which we will define in the following. Let (x,x r ) G R 6 . Let Xj and dx.j 
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be as defined after Q), and define 

Xj (^i, • • • , Xj —\, X ■ X j'+l 5 • • • 5 4 ?./V), 

(^l; ■ ■ • > %j — 1) ■£) + • • • ) —1; j *£fc+lj ■ ■ 

dkjfr = dxi ... dxj-idxj + i... dx k -idx k+l ... dxN- 
Then 7 ! and p 2 are defined by 

N 


x j = 


£ 


~ JW>N-3 


^(Xj)^(x') dxj, 


p 2 (x,x') = 

1 <j¥=b<N' 


R3N-6 


(x j:k )\ 2 dx jik . 


,xs). 


( 1 . 8 ) 

(1.9) 


In order to describe the regularity of 71 and p 2 we introduce D = 
{(ic, x) G M 6 } and 

S = ({R u .. ., R l } x M 3 ) U (M 3 x {R u .. ., R l }) C M 6 . 

Our method implies that 71 is real analytic on M 6 \ S and that p 2 is 
real analytic on M 6 \ (D U S). 


Remark 1.4. In the case of an atom, consider the density p in polar 
coordinates (x = ru>, r = \x\,u = x/\x\ G § 2 ), and define p{r) = 
J g2 p{ru)du. An important question is for which fcGN, 



( 1 . 10 ) 


exists. An even more demanding question is whether p(r) is real ana¬ 
lytic for r > 0 , i.e. whether p can be continued analytically beyond 0 . 
The analysis of such questions show the intimate relation between the 
problems i) and ii). In || it was shown that the derivatives in (|1.10|) 
exist for k < 2. But the general problem remains open. 


Remark 1.5 (Generalisations). As will be seen from the proof, The¬ 
orem 0 can easily be generalised to many other potentials. We do 
not use any of the special properties of Coulomb potentials, such as 
symmetry, homogeneity, etc. To be precise, let 

N N 

v ( x ) = & ) + S w jA x j~Xk), 

3=1 j,k=l,j^k 

satisfy the following conditions: 

(1) There exists C > 0 such that for all u G IT 1 , 2 (M 3JV ), 

II^ m IIl 2 (r 3jv ) < c , II m IIw i ’2(r3jv). 


( 1 . 11 ) 
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(2) There exists a constant L > 0 (depending on e) such that for 
all a G N 3 , we have 


N 


N 


5 ^ Ilcn/J + ]T Ww j:k \\ < i | “ |+1 | Q |!, (1.12) 

j = 1 j,k=l,j^k 

where the norms in (|1.12| ) are in L°° ({x G M 3 | \x\ > e}). 

The hrst condition, (|1.11|) , is a kind of relative boundedness assump¬ 
tion. The second condition, ( |1.12| ), means that V is real analytic away 
from £ (with a uniformity at infinity). Theorem |lTT| remains true for 
any V satisfying these two assumptions . For instance, replacing one or 
more of the Coulomb potentials in V by the Yukawa potential, e~ a ' x ' / |x| 
(with a > 0), we still get Theorem |1 . 1| . But here we concentrate on 
the physically important case of Coulomb potentials and do not strive 
for generality. 

Organisation of the paper: 

In section |3| we present in Lemma |3.1| a result concerning ‘partial an- 
alyticity’ of an eigenfunction if of H in the following sense: Upper 
bounds to the L 2 -norms of certain directional derivatives of if of arbi¬ 
trary order \a\ are given. They show the right behaviour in |a| needed 
later on for the proof of the analyticity of p away from the origin. (The 
proof of Lemma |3. 1| is given in Appendix |A].) We note that this kind of 
directional derivatives correspond, roughly speaking, to ‘taking deriva¬ 
tives along singularities of the potential’, see Lemma [A.3| and its proof, 
and compare also with PJ and 0]. Corollary 1333 is an immediate con¬ 


sequence of Lemma |3.1| and essential for the further steps in the proof 
of Theorem |1.1| . 

In section ft] we state and prove Proposition |4.1| which gives us the 
necessary control on |(<9 Q p)(;r)| for \x\ > e > 0 and arbitrary a. There¬ 
from the analyticity of p follows immediately. The key point of the 
proof of Proposition 0 is Lemma |4.3| . For its proof we use a suitable 
partition of unity of R 3JV and then proceed by a similar construction 
which together with Corollary 
Once Lemma 


as m 


ticular ( |4.10|) 
easy arguments. 


implies Lemma [L3 
is proved, Proposition 


m par- 
follows by 


2. Basic facts and notation 
Remember that for multiindices a = (oq,..., o 3 m) G N 3M , 

3 M 

M = ^2 a j- 

3 = 1 
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Furthermore, we have the usual ordering on multiindices: For a = 
(ai,..., o 3 m), (3 = b 3M ) we write a < (3 iff aj < bj for all j. 

We will need one simple and standard combinatorical fact. We recall 
it here for the reader’s convenience. 


Proposition 2.1. Let a G N 3JU be a multiindex. Then 


E 

P<a,m=b 



Proposition ^T] will be used as follows. Use Leibniz’ rule to calculate 

/3<a 


Then the number of terms where exactly b differentiations fall on / is 
given by ft 1 ). 

In the following we shall work with certain directional derivatives. 
Let e s for s G {1,2,3} denote the standard basis for M 3 . Let P be a 
(non-empty) subset of {1,..., N}. We define the coordinate xp by 


x P 


1 

TPl 


Ev- 

jeP 


We will now define d^ s p f for a function / G C 1 (M 3Ar ). For the given 
P and s let v = (v\,... ,Vn) G M. 3N with Vj = 0 for j ^ P, and 
Vj = e s j -ftP for j G P. Then we dehne 

9%. fix) = V/ ■ v. 

The dehnition of d® p then follows by iteration for any a G N 3 . One 
can clearly reformulate this dehnition in terms of Fourier transforms 
(multiplication by £J> S for suitably defined £p in Fourier space). In 
the previous paper [I]] we used a coordinate transformation to describe 
these derivatives. 


3. Partial analyticity of atomic eigenfunctions 

We will need a result on partial analyticity of the eigenfunctions of 

H. 

Lemma 3.1. Let -0 G L 2 (M 3JV ) be an eigenfunction of H . Let the 
index sets Pi,..., Pm C {1,..., N} satisfy for all s G {1,..., M } that 
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P s 7 ^ 0. Define for each s, Q s = {1, N}\P S . Define also, for e > 0, 


U P (£) = 


,..., Xn) G 


»3 N 


\xj | > £ for j G P s , 


\xj - x k \ > £ for j e P s ,k e Q s \. (3.1) 


Denote 


U Pl _ PM (£)=n*t 1 U Ps (£). (3.2) 

Then there exist C,L (depending on £) such that for all multiindices, 
a = {oil, ■ ■ ■ i °im) £ N 3M 7 we have 


llrf • ■■^J\\ mupi . Pm( «), < C£'“l( | Q | + 1) 


H 


The proof of Lemma [hi] is similar to the standard proof that solu¬ 
tions to elliptic equations with analytic coefficients are analytic (see || 
Section 7.5, pp. 177-180]) and will be given in Appendix |A|. 

Let us introduce the following practical notation. For a multiindex 
a = («i,..., %) € N 3M and given Pi,, Pm as in Lemma 
define <9“ p and t/p(e) by 


we 


d a = d ai 

% X Pl 


f)&M 


Up(e) - 


(3.3) 


We will need the result of Lemma 3.1 in a slightly different form for 


the proof of Theorem |1.1| . For later convenience, we state and prove 
this reformulation here. 


Corollary 3.2. Let the notation and assumptions be as in Lemma pJ 
(using or Then there exist constants Ci,Li such that 



rfrfwi At<crf' 

( £ ) 


(M + l) |a| 


Proof. By Leibniz’ rule, we have 
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Applying Cauchy-Schwarz and Lemma to both dxpi/j and <9" p 


m 


(|3~|), we find using Proposition |24] for the equality below, 




x) dx 


< cV w J2 (D (PI + i) M (H - PI + 1) W_W 

/3<a 'P' 

H /, IN 

= C 2 L^ J2 ( ? ) ( b + 1 ) 6 (l«l - b + 1) H " 6 

h —n x / 


6=0 

ll Q lf 


< C 2 (2L) H (|a| + 1)K 

Thus, Corollary fL2| holds with C\ = C 2 , L\ = 2 L. 


□ 


4. Differentiating the density 

Fix an arbitrary e > 0. We will always study p[x i) in the region 
{|xi| > e}. We will prove the following estimate: 

Proposition 4.1. Let e > 0 be given. Then there exist constants 
C, L > 0, such that for all |x| > e and all a E N 3 , p satisfies 

\d Q p{x)\ <CL H (|a| + l) H . (4.1) 

Remark 4.2. It is clear that Proposition |4.1| implies Theorem |1.1| . 

Proof of Proposition m Choose xi,x 2 £ C 00 ^ 3 ), satisfying 

Xi + X ‘2 = 1, XI = 1 on #(0, e /( 4Ar )) ; SU PP Xi c B( 0, e/(2N)), 

and let further XhX 2 be radially symmetric functions. Using this par¬ 
tition of unity and the notation 




j <k\ , 


we can write 
p( x i) = 


x )| 2 n (Xl ( x 3 - X k ) + X2 {Xj - x k )) dx2 • • • dx N 

j<k 

/ L /. 12 / 


Y. / |V’| 2 (x)0/(x) dx 2 ■ ■ -dx N 
TcM J 

Y pi( x i) 


(4.2) 


ICM 
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Equation (f4.2|) defines 0/ as 
4>i = 


n 

U,k)ei 


Xj 


n 

(j,k)eM\I 


iCfc)}- (4-3) 


We will prove that pi(x \) satisfies an estimate like (|47T|) on {|xi| > e} 
for all / C Af, namely 


\d a Pl \ <CL^ (|a| + l) |a| . 


(4.4) 


The estimate (|4.1|) follows from (|4.4|) (with a different C) since the 
sum in ( |4.2| ) is finite. 

The estimate ( |4.4p is a consequence of ( |4.1(J| ) (with 0 = 0/) in Lemma 
|4.3| below, using a Sobolev imbedding theorem. Since we have not found 
an ideal reference we include the following easy argument: 

Let v G (^“(M 3 ), v(x) = 1 for |x| > e, v(x) = 0 for \x\ < s/2. Let 
furthermore T denote the Fourier transformation. We get for a G N 3 , 
|x| > e, 


d a p(x) = (vd a p)(x) 

= cj e ixp {l+p 2 )~ 2 P ((1 - A) 2 (■ vd a p )) dp. 

Therefore, 

\d a p(x)\ < c||(l +p 2 )- 2 || il(R 3 ) ||(l - A) 2 (u9“p)|| L i (R 3 ) . (4.5) 

We can estimate ||(1 — A) 2 (ucP//)||li( R 3) using ( |4.1Q ) by, 

11(1 - A) 2 ( ! ,a»|| iI(E , ) < Cl 4“ l+4 (|a| + 4 + l)'“l+ 4 

+ (4.6) 


for some constants ci, L\, C 2 , L 2 . Combining (f4.5|) and (|4.6|) yields ( |4.4|) . 
Proving Lemma |4.3| therefore finishes the proof of Proposition |4. 1|. □ 


Lemma 4.3. Let e > 0 be given and let 

4> = II fj,k( x j~ x k ), (4.7) 

1 <j<k<N 

where each fj k is one of the functions XhX 2 ,d es X 2 , with e s G N 3 , 

|e s | = 1. 

i) Let Pi, ..., Pm be subsets of {1,, N } satisfying that 1 G Pj 
for j = 1,..., M and 

(supp^) n {|ii| > e} c n" 1 (7p f (e/(4Af)). 


(4.8) 
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Then there exist constants C,L > 0 (depending on e) such that 
for all multiindices, a = (a 0 , ol\, ..., ckm) G N 3M+3 , we have 



■ ■ ■ d xp M \M 2 ) ( x )0( x ) dx 2 • • • dx N 




<CL h (H + 1) |q| . 


(4.9) 


ii) There exist constants C, L > 0 (depending on e) such that for 
all a G N 3 we have 


d a 

X\ 


(x)0(x) dx 2 ■ • ■ dx jv 


^({kil>4) 


<CXl“l(|a| + l)K (4.10) 


Proof. To a function f> as given in ( |4.7|) we will associate P = P(f>) C 
{1,..., iV} satisfying 1 G P and such that 


(supp0) n{|xi| > e) C t/p ( 0)(e/(4iV)). 


(4.11) 


We will now describe the map (j) i—> P(<f>). We note that the following 
construction is similar to the one from [jlj. Define / = 1(f)) C by 

(j, k) e I(f>) if and only if f jjk e {xi, d ei X 2 , 9 e2 y 2 , 9 e3 x 2 }- 

In other words, (j,k) G 1(f)) means precisely that fj ik ^ y 2 . The set 
I(f>) generates an equivalence relation on {1,...,1V} 2 and we define 
P(f>) to be the equivalence class of 1. Less abstractly, this means that 

• 1 e P(f)). 

• For j > 2 we have j G P(f>) iff there exists {.)i , • • •, j.s } Q 
{1 .....iv}. s < N, satisfying 

(l)Ji) ^ d(ft), 

(jt,jt+ 1 ) G 1(f)) or (j t+ i,jt) G 1(f)) for 1 < t < s, 

and 

(js,j) e 1(f)) or (j,j s ) G 1(f)). 

Notice that, since yi+x 2 = 1, snpp d ej X 2 C snpp xi, j = 1, 2, 3. There¬ 
fore we get ( |4.11|) by the same elementary geometrical considerations 
(the triangle inequality) as in [1J. 

In the proof we shall use P = P(f>) in order to replace the derivative 
df° outside the integral in the left hand side of (fO| ) by the derivative 
df° p inside the integral. That will enable us to apply Corollary |3.2| . 

Let P = P(f>) according to onr construction. We will prove the 
lemma recursively in |P|. In the proof below we will freely interchange 
the order of differentiation (in the distributional sense) and integration. 
This is permitted, due to Corollary |3.2| , which ensures that the deriva¬ 
tives of the functions in question belong to L 1 ({|aq| > e} x R 3Ar-3 ). 
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We will only prove part i) of Lemma |4.3| . The changes necessary for 
the case ii) are obvious and therefore omitted. 

Step 1, \P\ = N. In the case where P = {1,...,1V} we make 
the change of variables i/j = Xj — Xi for j = 2, ..., N. Then we get 
Xj — Xk = Vj — Vk for j, k 7 ^ 1. The point is that 0 only depends on 
the differences Xj — Xk, and therefore, after the change of variables, the 
only dependence on x\ will be in |0| 2 , where we can apply Corollary 
p. 2 | . Let us carry this out. 

Denote y = (y 2 ,..., Vn)- Then we see that after change of variables 
we have 

0(x) = 0(y), 


for some function 0. Explicitly, we see from (14.71) that 


N 


0(2/2, ■ ■ ■, vn) = n fjAvj - n AA-Vs)- 

1 <j<k 


s=2 


Therefore, 




|Ql 


l'0| 2 )( x )0(x) dx 2 ■ ■ ■ dx N 


(4.12) 


= / ( dxl , ■■■d^J^\ 2 )(x 1 ,x 1 +y2,...,x 1 +y N )(t>(y) d y- 


From ([4.12|) we get by differentiation under the integral sign and change 
of coordinates back to x: 


/ (dT Pl ■ ■ ■ d xp M l^| 2 )( x )0( x ) dx 2 • • • dx N 


— (pjOLQ p)OL 1 ^ # . pjOLM 

/ \ U xp U xp 1 X P M 


(4.13) 

\ 2 )(xi,x t + y 2 ,...,x 1 + yn)4’(y)dy 
= J (BZSZl, IV>| 2 )(x)^>(x) dx 2 ■ ■ ■ dx N . 

Notice the support conditions (|4.11|) , ( |4.8|) . We can now apply Corol¬ 


lary |3]2 to get ( [h9D in the case \P\ = N. 

Step 2, | P\ < N. Suppose that Lemma |4.3| holds under the addi¬ 
tional assumption \P\ > K for some 0 < K < N. We will prove the 
statement for \P\ = Ii. 

Dehne Q = (1,..., IV} \ P. Since |P| < N, Q ^ 0. Note that if 
j e P, k G Q then, by definition of J(0) and P = P(0), we have 
(j, k) 0 I (0) and ( k,j ) 0 7(0)- Therefore, if j < k we have fj t k = X '2 
and if k < j we have fkj — X r 2 - So 0 contains the factor (remember 
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that \2 is rotationally symmetric, in particular even) 


b,Q = JJ X 2 (Xj-Xk), 
jeP,keQ 


and can be written as 


where 


Op ■ 


0P | I j, kip j %k)i 0Q | fj,k (pj 3Ck)- 

j,k£P,j<k j,k£Q,j<k 

We do the following change of variables for j > 2: 

J Xj - Xi for j E P \ {1}, 
yi= \ Xi for j € Q. 

For convenience of notation we define y\ = 0. We clearly get for j,k E P 
or j, k E Q, that x 3 — = Vj — ])k (also when either j — 1 or k — 1— 

remember that 1 G P). Thus, as in the case \P\ = N, we can write 
<M X ) = 0p(y), and <fo(x) = <fo(y). 

j = 1, 

Write z = (^, z 2 ,..., z N ) E M 3W with ^ xi + yj j E P\ {!}, 


Xi 


J e Q- 


Then 


(&xL ■ ■ ■ d xi* H 2 )( X )<K X ) 


(4.14) 


= / •••^M 2 )( z )<My) II wOd + %- -yk )) 4>qi.y)dy. 

J \jeP,k£Q 

Differentiation under the integral sign yields 


V,, / (^ • • • d™ M |^| 2 )(x)0(x) dx 2 • • • da* 


(4.15) 


= / (V, P ^---^|^)(z)0p(y){n wOd + yj -yk)}4> Q (y)dy 

d jeP,keQ 

+ E [Pl,--P?JP)(*)My)My) 

j£P,kGQ ' 

x( n X2(®i + J/j'-J/fc'))(Vx2)(a:i + J/j-J/fc)dy. 

j'eP,k'eQ,{j',k')^{j,k) 

Let us explain roughly how we will proceed for higher derivatives with 
respect to X\. For each consecutive differentiation we will get terms as 
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in ( |4.15|) . The term where all differentiations fall on t/j can be differen¬ 
tiated again in a manner similar to ( |4.15| ). If one differentiation falls on 
X 2 we stop differentiating that term under the integral sign—leaving 
the rest of the differentiations outside the integral. The result of this 
procedure is (|4.16|) , the notation of which we will define below. The 
important point is that when all derivatives fall on i/j, we can apply 
Corollary fl.2| to obtain our conclusion. On the other hand a differen¬ 
tiation of X 2 will lead to a situation with a larger \P\ —so these terms 
can be handled by the induction hypothesis. 

Let 7i s G N 3 , \r} s \ = 1, s = 1, 2,..., S. Define for t G {1,..., S} 

s 

4) = 0 , A = 5^’ Bt = ^2 Vs, B s = 0 . 

s<t s=t -\-1 

Notice that the definition of B t depends on S, i.e. B t = B t (S). We get 
the following formula ( [4 .1 6| ) from (|4.15| ), using the procedure described 
above, by induction with respect to S. For 5 = 1 the equation ( [4 .1 6| ) 
reduces to (|4.15|). 




(4.16) 


J {dxp d T Pl ' ' ' 8 xp m l^l 2 ) ( X )<K X ) dx * ■ ■ ■ dx N 

s r 

+ E 9 ®.‘{ E y)*j(y) 

t =1 ieP.keQ J 


j£P,k£Q 
X 


n 

j'eP,k>eQ,(f,k')^(j,k) 


A 2 A 1 + yy ~y k ')) 


x {d Vt X 2 )(xi + yj - y k ) d V 2 ■ ■ 

We will use ( |4.16|) with S = |a 0 |j As = a 0 . Consider the function 
<l>j,k = <i>P-<i>Q-(d nt X 2){Xj-X k ) JJ X2{Xj'-X V ). 

By construction we have \P((f>j,k)\ > |-P(0)|- Therefore, we get by the 
induction hypothesis on \P\ that 




J • • • d xp M IVf)( x )<M X ) dx 2 ■ ■ ■ dXN 

<CL p { P +iy, 

where p = \B t \ + |A-i| + |«i| + ... + \oim\ = M — 1. Furthermore, 
using Corollary |3. 2| on the first term on the right hand side in (|4. 16|), 
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we obtain 



'' • d xp M |-0| 2 )(x)0(x) dx 2 ■ ■ ■ dx N 


< cisip+iy, 

LH{\xi\>e}) 


with p = |Ag| + |a;i| + ... + |CKjvf| = |a|. 

Thus the desired estimate holds for the individual terms on the right 
hand side in ( |4.16|) . Since the number of terms is bounded by c|a| this 
finishes the proof of (|4.9| ). □ 


Appendix A. Proof of Lemma 3.f 


In this appendix we will prove Lemma fi. 1| . For convenience define 
He = H — E, with E being the eigenvalue corresponding to the eigen¬ 
function 0, i.e. 0 satisfies He 0 = 0. Recall the notations given in 

(ED- 

Let us start the proof by stating a well known result explicitly. Since 
the domain of He is known to be IF 2,2 (M 3Ar ), we get 


Lemma A.l. Let v G kF 1,2 (M 37V ). Then v G IF 2,2 (R 3Af ) if and only 
if H e v G L 2 (M 3Ar ). Furthermore, there exists a constant K 0 > 0 such 
that for all v G PF 2 ’ 2 (M 3JV ) 

|M|w 2 ’ 2 (R 3JV ) A A 0 (\\H E v\\ L 2 {R 3N ) + ||u||i2( R 3JV)) . 


This follows from the fact that V is infinitesimally small (in the 
operator sense) with respect to —A. 

We now state and prove an a priori estimate. 


Lemma A.2 (A priori estimate). There exists a constant Co such that 
for all jp jj i G (0,1), all a 0 G N 3N , a' G N 3M with |a 0 | + |a'| < 2 and 
all v G W 2 ' 2 (U F (£/2 + rh)) we have the estimate 

>;l“ol +l “'l||9"»S I ">l|y ([ / p( « /2+1J+ , l)) (A.l) 

<C'o{l? 2 ||ffEl>||Ll([/,(«/2+, 1 )) + E, ’l l ' 3| |l a ' 3t ’||i>([/p( e / 2+ „))}' 

/3eN 3Ar ,| / g|<2 


Furthermore, if the right hand side of (]A . 1|) is finite for all rj, rp > 0 
then v G W 2,2 (Up(e/2 + rft)) for all p\ > 0. 


Proof. Since Up (e/2 +rj + rji) C Up(e/2 + pi), the estimate is obviously 
true for |a 0 | + \ a '\ < 2. Let a 0 G N 3Ar , of G N 3M with |a 0 | + |cd| = 2. 
Choose <fr G C yoo (M 3Ar ), 0 < 0 < 1, with 0=1 on Up(e/2 + rj + p x ) 
and supp0 C Up(e/2 + rji), satisfying ||0 7 0|| oo < with C 7 

independent of rj,rji. 
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We can now estimate, using Lemma in the third inequality below 
\\ dfX ° d x p v\\ L 2 ( u P (e/2+ V + Vl ) ) < ||5 a °^p(HllL2 ( R3 J V) 

< c||(0u)|| W -2,2( R 3iV) < C (||i/£;(0T)|| L 2( R 3JV) + ||0U|| L 2( R 3JV)) 

< c|||0i/£;u|| i 2( R 3 iV) + ||(A0 )u|| L 2( R 3JV) 

+ 2|| V0 • Vull^RSiV) + ||||^2( R 3Af) j’ 


< C 


•P/ 2 + 771 )) 


+ Ci?7 \\Vv\\ L 2^j p ^/ 2+m ))+c 2 r) \\v\\l2(Up{£/ 2+ Vi ))\, 


for some constants c, C, Ci,c 2 . Inequality O follows by multiplying 
with if. 

The last statement of the lemma follows easily from Lemma |A.1|. □ 


Finally, we state and prove the properties of V that we need in the 
proof of Lemma |3.1| . 

Lemma A.3 (Properties of V). Let V be the Coulomb potential defined 
in (|L7D. 

(1) There exists C v > 0 such that for all v G kF 1,2 (M 3iV ) we have 

II iy — -^)t||l 2 (r 3JV ) < Cy ||n|| w 1 - 2 (r3jv). (A.2) 

(2) There exists a constant L v > 0 (depending on e) such that for 
all a G N 3M with |a| > 1, we have 

ll^pr|L» ([ /p(«/2))<il? l+, |Q|!. (a. 3 ) 

Remark A.4. These are the only properties of V that we need (to¬ 
gether with Lemmas M and |A.2| ). They are easily seen to hold (see 
argument in proof below) for potentials satisfying the general condi¬ 
tions in Remark |1.5| . 


Proof. The first property (|A.2|) is a consequence of Hardy’s inequality 
(see for instance || Vol. II, p. 169]). To prove the second property, 
(|A.3|) , let P s be one of the index sets defined in Lemma |3.1| . Notice 
that 

0 for j f P s 


| xA~ l = 


and 


'x Pl 


|P,l“ |a| / 2 a?b| _1 


for j G P s , 


dfpJxj-Xk | = 


0 

|p |-H/ 2 d°:\x\~ l 


X=Xj—Xk 


for j, k£P s or j, k G P s 
for j G P s: k (f P s . 
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Therefore, ( A.3 ) follows from the structure of V, the real analyticity of 
x i—> |x| _1 away from 0 and the definitions of Up {s/2) and df . □ 


Proof of Lemma \ 3. 1\ . Notice that ( |A.3|) trivially implies that for j, r] > 
0, jrj < 1, |a| > 1 we have 




(A.4) 


We will prove that there exists > 0, such that for all r] G (0,1) 
and all j G N with jtj < 1 we have, for all a G N 3M , «o G N 3JV , «o < 2, 
| ot | + |do | < 2 + j, 

,H+KIn^o< L H+I“»I«. (A.5) 

Before proving ( |A.5| ), let us note that Lemma |3j] follows easily from 
it. In fact, let a G N 3M , | a\ > 1 and choose |a 0 | = 0, r) = e/(2|a|), 
j = a|. Then ( |A.5|) becomes 

l|3“p^lli>( [ /p(A)<h“ l+1 (2/e) W |a| l “ l . 


which is the statement of Lemma 3.1 


We now prove (|A.5| ) by induction in j. For j — 0,1, there is nothing 
to prove since we know that if G IF 2,2 (R 3Ar ). Let L^ be sufficiently 
large for (|A.5|) to be true for j = 0,1 and satisfying furthermore, 


> max 12 L v , C 0 (l + ^ : 

\3\<2 


(A.6) 


Here the sum is over all Q G N 3JV with \f3\ < 2, Cq is the constant from 
Lemma [A.2| , and Ly is the constant from ( |A.3| ). 

Suppose that we have proved (|A.5|) for all j < j 0 and all rj G (0,1) 
with ijj < 1. We will prove that (|A.5|) holds for j = jo + 1 and all 
7] G (0,1) with (jo + 1)77 < 1. 

Let |a| + |a 0 | < 2+j 0 . Then clearly Up(e/2 + (j 0 + l)r]) C f/p(e/2+j 0 ?/). 
Therefore, 


v' al+lao '\\d ao d ^\\ L 2 (UT>(e/2+( jo +i) v )) <V lal+laol \\d ao df p 7f\\ L 2 {Up(£/2+jori)) , 


and the result holds by the induction hypothesis. So we only have 
to consider the case |a| + |a 0 | = 2 + j Q . Choose a decomposition 
a = of + a", with \a!\ = 2 — |a 0 |, be. with \a"\ = j 0 . Using Lemma 
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A.2| with r/i = j 0 r), v = <9" p ip we find 


7/ 2+Jo 11 8 ao <9 “ p Ip 11 L 2 (c/p (e/2+Oo+l)^)) 

< C 0 \r) 2+jo \\H E d: 


\ L 2 (Up (e/2+jorfj) 


E 7 7 l/3|+io 11 ^ I1L 2 (Cp (e/2+ion)) } ■ ( A -7) 


|<2 


Since H E = 0, we get 

r/ 2+i0 11 tf E d“ p V|| L 2 ([/ P (e/2+ior;)) 


= 7^~b70 


E 

7<ct",/3+7=a" 


a 


7 


(9f p (y-£)K> 


L^U F (e/2+j 0 rj)) 


(A.8) 


We now use (|A.4f) , the combinatorical result from Proposition ^]l] to¬ 
gether with the induction hypothesis, to estimate ( A.8 ) as 

V 2+jo \\H E d^\\ L 2 {Up(E/2+jov)) 

KI 


< 


Kl /, //, 

a 


k 


E 

fc=i 


< Y. L v +lL 


/c-f-1 T D7 - ^+1 


k =1 


< <’ l+2 < hf l+2 - 


(A.9) 


k =1 


Here we used the assumption > 2 Ly from ( A.6 ) in the last estimate. 

Due to the induction hypothesis we can also estimate the other term 
in (|OD 

E ^ii^Oii< E L *" I+W+1 - < a - 10 ) 

l/3|<2 |/9|<2 

So using (|A.9| ) and ([A. 10|) , we can estimate ( |A.7p as 

2+jo II f) a ° f) a i/i|| „ <7 rM + l“o|+l / A ^-d/ 3|<2 ) 

V .x P Yl|L 2 (I/ P (e/2+0'o+l)D) — 


The last factor is < 1, by the choice of (see (|A.6|) ). That finishes 
the proof of (|A.5|) and therefore of Lemma |3.1|. □ 
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